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Abstract. 

Compelling evidence exists to show that the structure of molecular clouds is 
fractal in nature. In this poster, the author reiterates that view and, in addi- 
tion, asserts that not only is cloud geometry fractal, but that they also have a 
common characteristic - they are similar in shape to the Horsehead in Orion. 
This shape can be described by the Julia function f(z) = + c, where both z 
and c are complex quantities and c — -0.745429 4- 0.1130081. The dynamical 
processes responsible for the production of these clouds seem to be turbulence 
followed by Brownian motion till high densities arc reached, at which point 
structure formation is dictated by gravity. The author presents image analysis 
of four varied examples, namely those of the Horsehead nebula. Helix neb- 
ula, Eagle nebula. Rosette nebula and Paley I nebula to prove her hypothesis. 
The images of fractal these nebulae are analyzed for their box dimension using 
fractal analysis software and comparisons are made with the given Julia set. 
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1. Introduction 

Molecular clouds are self-gravitating, magnetized, turbulent, compressible flu- 
ids (J.P.Wilhams et.al. 1999; Elmegreen 1999; R.B.Larson 1995). In addition, 
their molecular structure is similar over a wide range of scales with similar 
power law indices, quite independant of the star-forming nature of the cloud, 
even at a resolution of a few parsecs. The only difference between these clouds 
has been found to be in their average densities in that in regions of high den- 
sities, gravity begins to dominate and self-similarity breaks down, ultimately 
leading to the formation of stars, planets, etc. This implies a hierarchy of 
structure, beginning from the giant molecular clouds with sizes about 50 pc 
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and masses ranging from 10"* to 10^ AfQ.Thc next class is the normal clouds 
with sizes ranging from a pc to 30 pc and masses from 25 Mq to 10^ Mq 
(S.S.Prasad 1995). These clouds may or may not be bound and may contain a 
small number of low mass stars. However, their contribution to the total star 
formation rate in the Galaxy is negligible. The third is the clump within the 
cloud. Clumps have sizes ranging from 0. 4 to a few pc amd masses from a few 
to 25 Mq. They are coherent star-forming regions out of which stellar clusters 
form. The gas clumps are normally bound by pressure rather than by gravity. 
However, both clumps bound by strong self-gravity and by external pressure 
have identical density distributions. The fourth and final class is that of cores 
found within clumps, which have sizes 0.1 pc tp 0.4 pc and masses less than 
1 Mq. Cores are regions out of which single stars (or multiple systems) are 
born and are always gravitationally bound. Stars formed in cores may accrete 
matter from the host clump or cloud at the protostar stage. 

Power-law relationships in molecular clouds have been found in earlier quanti- 
tative studies (Kramer ct. al.l998; Hctcm & Lcpinc 1993;Stutzki ct.al.l998). In 
this poster, the author investigates the structure of certain molecular cloud im- 
ages, by using a fractal dimension analysis software. The images presented are 
of (i) Horschcad nebula (B33) in Orion (Kramer ct. al. 1996), (ii)Eaglc nebula 
(M16) in Serpens (Hester et.al. 1996), (iii)Rosette nebula (NGC 2237 to NGC 
2246) (Carlqvist et.al. 1998;Clayton et.al. 1998; White 1997) and (iv) IRAS 
02356-2959 (Palcy I) in Fornax (Palcy et.al. 1991; Stark 1995). Of these four, 
that of the Horsehead nebula has been obtained from the author's observations 
from the 2.3m Vainu Bappu Telescope, Kavalur, India. Basic mathematical 
backgroimd for this work is presented in section 2. , observations in section 3., 
analysis in section 4. Finally, discussions are presented in section 5. 

2. Mathematical Background 

Fractals are self-similar objects, having fine structure. Examples of fractal ob- 
jects include the von Koch curve, the Sierpinski gasket, the Cantor set and 
the Julia sct( K. Falconer 1997; H.O. Peitgcn et.al. 1992). In order to study 
fractal structure in physical objects, we need a quantitative measure defined on 
them. Such a quantity is the principal definition of fractal dimension , other 
definitions, such as Kolmogorov (or Minkowski) dimension and box (or grid) 
dimension, are widely used. Indeed, in practise, box dimension is the easiest to 
work with and is equivlent to Kolmogorov dimension. The box-counting dimen- 
sion is defined in the following way : if X is any non-empty bounded subset of 
R" and N^{X) be the smallest number of sets of diameter at most e which can 
cover X, then box dimension is defined as 

dimX = \im logN,(X)/log(l/e) 

In the above definition of dimension, the fundamental assumption is that for 
each e, we measure a set in such a manner that irregularities of size less than 
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e will be ignored. For example, if X is a plane curve, then Nf^{X) will be the 
number of steps required by a pair of dividers set at length e to traverse X. We 
have to also keep in mind that for real phenomenon we can use only a finite 
range of e. An example of a fractal object is that of the Julia set given by the 
function 

f{z) = z^ + c 

where z and c are complex numbers and c = -0. 745429+0. 113008i(Fig.l) How- 
ever, unlike Euclidean objects, a fractal object retains a similar level of com- 
plexity for each value of e, the fractal dimension gives an index of the degree to 
which a fractal structure fills the space in which it is embedded (B.Mandelbrot 
1982) Hence, examination of an object for self-similarity at different scales of 
magnification can suggest whether the object has a fractal structure. 

3. Observations 

Images of the Horsehead nebula were obtained with the 2.3 m Cassegranian 
telescope at the Vainu Bappu Observatory (VBO) located in Kavalur, India, 
1*** and 3'''^ April 2000. The SITe 1024x1024 pixels, 24 micron charge- coupled 
device (CCD) was placed at it's prime focus. The camera projects each pixel 
to 0.66" on the sky (resolution of 27.6" per mm) giving a field of view of 11.264 
X 11.264 ". The images were obtained through H-alpha (6563 A)(Fig.2), V 
and R filter with FWHM of H-alpha filter being 100 A, while the V and R 
filters were Bessel filters with the standard bandpass. The effective exposures 
of the images are given in Table 1. The seeing during the observation of the 
order of a few arc seconds. Each image was bias subtracted and flat-fielded 
using twilight fiats obtained each day. However, dark currents were assumed 
to be neglible as the camera was cooled to about - 90° C. Both bad regions 
and cosmic rays were removed using IRAF NOAO image reduction software. 
The CCD, with a read-out noise of the order of 8 electrons, was calibrated 
using IRAF software. The standard stars used for guiding the telescope were 
obtained from the Hubble Space Telescope (HST) Guide Star Catalogue. 



Table 1. Observation details from the 2.3m VBT, Kavalur. 
Co-ordinates are based on year 2000. 



Date 


Right Acsension 


Declination 


Filter 


UT 


Exposure time 


2.4.00 


5h40m58.6s 


-2'J27'23.84" 


H-alpha 


1357 


30 


3.4.00 


5h40m58.6s 


-2"27'23.84" 


R 


1353 


15 


3.4.00 


5h40m58.6s 


-2"27'23.84" 


V 


1426 


30 
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4. Fractal Analysis 

The images (Figures 3 to 5 ) were obtained from various sources (HST web- 
site; D.Malin 1983;R. Stark 1995) and converted to Windows bitmaps. This 
causes the digitized image to be stored as a data matrix where pixels be- 
longing to the pattern are stored as 1 and pixels from the background are 
stored as 0, or vice versa. This procedure can also be carried out for colour 
images (which are overlays of 3 grey-scale images. The fractal dimension of 
these images were then measured using an implementation of the box-counting 
method. The automated fractal image analysis software, Benoit 1.3( Truesoft 
International Inc., St. Petersburg, USA), was used to convert the solid images 
to black and white pixel outlines. This software has been reviewed( Science 
1999, 285, pg.l228)and has been found to perform satisfactorily. These images 
were overlayed with grids in such a way that the minimum number of boxes 
were occupied. The is accomplished in Benoit by rotating the grid for each size 
through 90 degrees at various angular increments of rotation. The user is able 
to change the box length decrease factor as well as the size of the largest box. 



Figurel.The Julia set for c = -0.745429-1- 0.113008i(Box dimension= 1.679594) 



For the analysis, a set of ten measurements were taken at fixed grid rotation 
increments of 15 degrees. An average dimension was calculated and the number 
of iterations giving the slope with the least standard deviation (SD) while being 
nearest to the to the average was chosen. Log-log graphs were then plotted of 
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Figure2. H-alpha image of the Horsehead nebula obtained with the 2.3m VBT, Kavalur(Box 
dimension=1.6965725) 



Figure 3. Part of Paleyl (Stark 1995) (Box dimension= 1.67718) 



Figure 4. Closeup of a trunk of the Rosette nebula (Malin 1983)(Box dimension=l. 78261). 



the reciprocal of the side length of the square against the number of outline 
containing squares. The slope is estimated by fitting a line using the method 
of least squares. Ten further sets of measurements were taken with this selected 
number of iterations and box size decrease factor. Results for each cloud were 
tabulated and then the Student's t test of significance (Table 2) was used 
to test for their deviation from the value of the dimension for a Euclidean 
shape ( topological dimension 1). A further test was made for fluctuations of 
the cloud dimension values from that of the Julia set. The results show that the 
cloud dimensions are significantly different from it's topological dimension at 1 
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Figure 5. Closcup of a trunk of M16(part of column I in Hester 1995)obtained from the 
Hubble Picture Gallery(dimension=l. 669622) 



percent cofidence limit (t value 3.25) but not from the Kolmogorov dimension 
of the Julia set. Hence, it can be said with certainity that molecular clouds are 
not only fractal but also that their dimension is that of the given Julia set . 



Table 2. Analysis of box dimension . 



Name 


Pattern size(pixels) 


Dimension 


Student's t value | 


Julia set 


794x1324 


1.679594 


91.19396 


Horsehead nebula 


406x473 


1.6965725 


16.52918 


Palcy I 


463x442 


1.677181 


95.62187 


Rosette nebula(close-up) 


109x81 


1.78261 


85.69697 


Eagle nebula(close-up) 


188x212 


1.669622 


70.454 



5. Discussion 

The discovery of power-law relationships in molecular clouds by Larson(1981) 

was a ground- breaking study in favour of self-similarity, especially so since a 
power-law can only be applied to clouds with masses greater than 10'' Mq and 
does not apply to cores and gravitationally bound clumps. My study on the 
different types of nebulae reinforces this view, especially since the parts studied 
are small sections of the main trunk in some clouds (eg. Eagle and Rosette). 
What is significant here is that their dimension is the dimension of the Julia 
set. This result implies that physical features of clouds can be explained if 
we make the following assumptions (i) cloud structure resembles the Julia set 
corresponding to the Mandelbrot constant given by c = -0.745429-1-0. 1130081 
and that it is a dynamical entity composed of cloud matter, (ii) growth occurs 
around the central point (an attractor) of each spiral making it a potential 
site of star formation, (iii) groups of such central points spiral inwards towards 
a massive central point, identifying such central points to be dark objects. 
Such central objects may become supermassive black holes (Macchetto et. al. 
1997) or remain as numerous small objects( Goodman and Lee 1989), perhaps 
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even composed of strange matter (van der Marel 1997; Weber 2000; Chandra 
& Goyal 2000) (iv) central points increase in size as we move away from the 
massive central attractor, (v) growth also occurs along each spiral of the set 
in the form of horseheads, whose sizes reduce adinfinitum as the spiral moves 
towards the attractive fixed point and finally (iv) after the elapse of sufficient 
time, when adequate matter has accumulated in the centre, a dwarf, a star or 
group of stars would form, depending on the size of the clump. 
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